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We have developed a two-electron outer region for use within R-matrix theory to describe dou¬ 
ble ionisation processes. The capability of this method is demonstrated for single-photon double 
ionisation of He in the photon energy region between 80 eV to 180 eV. The cross sections are in 
agreement with established data. The extended RMT method also provides information on higher- 
order processes, as demonstrated by the identification of signatures for sequential double ionisation 
processes involving an intermediate He^ state with n = 2. 


The development of laser sources capable of generating 
ultra-short light pulses or intense XUV and X-ray laser 
light m has put additional emphasis on the investiga¬ 
tion of multiple ionisation processes [4]. To complement 
recent experimental developments, there is a need for new 
computational techniques able to describe the multiple 
ionisation of general atoms. Although progress has been 
made in the description of double photoionisation of gen¬ 
eral atoms [S] , further progress is needed to enable inves¬ 
tigation of these processes in a time-dependent manner. 
This is of particular importance in experiments involving 
ultra-short light pulses. 

One possible route to develop capability for the de¬ 
scription of double-ionisation processes is R-Matrix the¬ 
ory. Time-independent R-matrix theory has already been 
extended for this purpose. The R-matrix with pseudo¬ 
states technique (RMPS) utilises a large basis set in the 
inner region, which includes residual-ion states in the 
continuum. Although the first electron can escape, the 
second must remain bound. Nevertheless, information on 
double ionisation can be obtained through evaluating the 
probability that the residual ion is left in an excited state 
above the ionisation threshold [6]. A second R-matrix 
approach for double ionisation is intermediate-energy R- 
matrix (lERM) theory, in which two electrons are allowed 
to escape the inner region. This leads to a 2-dimensional 
propagation of the R-matrix over a large distance where 
the R-matrix is matched to asymptotic solutions. This 
approach was initially applied to scattering 13IH], and 
has recently been applied to double photoionisation [^. 

One promising computational method for solving the 
time-dependent Schrodinger equation (TDSE) for many 
electron atoms in laser fields is R-Matrix with time- 
dependence (RMT) [TUHH]. With electron-electron in¬ 
teractions fully incorporated, RMT has been used to suc¬ 
cessfully model a number of features of many-electron 
atoms undergoing photo-ionisation Unilll]. In RMT, an 
R-Matrix basis set inner region m is attached to an 
outer-region finite difference grid [16] with the ability to 
describe an ionised electron, enabling the physical pro¬ 
cesses in each region to be modelled with an appropriate 
numerical method. 

While previous applications of RMT to double¬ 
ionisation m have calculated accurate two-photon cross 
sections, the range of observables obtainable with this 


method is narrowed by an upper limit on the radial dis¬ 
tance travelled by the inner ionising electron, causing re¬ 
flections upon interaction with the inner-region bound¬ 
ary. In this communication, we report a method capa¬ 
ble of modelling both ejected electrons accurately over a 
broad range of energies and integration volumes. 

In the present approach, we adopt the philosophy of 
lERM theory within RMT theory by allowing two elec¬ 
trons to escape into the outer region. Electron i in the 
outer region is here described through a finite-difference 
(ED) representation for the radial coordinate r^. We 
thus employ three distinct regions: (1) an inner region, 
in which all N electrons are within a distance a from 
the nucleus, represented with a standard R-matrix basis 
set, (2) a one-electron outer region, in which we com¬ 
bine a basis set representation for the residual system 
with {N — 1) electrons, and a finite-difference represen¬ 
tation for the ionised electron with > a, and (3) a 
two-electron outer region, in which we use a basis set 
representation for the residual system with {N — 2) elec¬ 
trons and a finite-difference representation for the two 
electrons with rpf-i,rN > a- 

As a proof-of-principle, we apply it in the present com¬ 
munication to single-photon double photoionization of 
He. This process has already been studied extensively, 
both experimentally m and theoretically 53110 ], and 
it is therefore possible to compare the outcomes of the 
calculations with benchmark data. However, since the 
approach described here allows the simultaneous calcu¬ 
lation of higher-order processes (such as double electron 
above threshold ionisation [2T|), we also expect signa¬ 
tures of these higher-order processes to appear in the 
results. 

The Hamiltonian for the helium atom in a laser field 
is given in atomic units as 

H = - 1 - \-E{t)-{zi+Z2) ( 1 ) 

2 2 ri r2 ri2 

and the time-dependent Schrodinger equation is given by 


In these equations, ri and r 2 are the radial coordinates 
of the first and second electron, and ^ represents the 
inter electron repulsion. Eit) is the time-dependent laser 
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field, and is the two-electron wavefunction. zi and 
are the positions of the electron in the direction of the 
laser field. 

The two-electron outer region wavefunction is de¬ 
scribed on a set of two-dimensional FD grids given as 
*" 2 ) where each grid describes the component of 
the wavefunction corresponding to the angular momen¬ 
tum quantum numbers across the coordi¬ 

nates ri = a...b and r 2 = a...b. Near the boundary with 
the one-electron outer region (ri « a and r 2 ~ a), the 
two-electron outer region is provided with the necessary 
wavefunction information from the one-electron outer re¬ 
gion on an extension of the FD grid. This allows the 
direct evaluation of equation ([^ using FD techniques. 

In the one-electron outer region, the inner electron is 
described using a near-complete basis of eigenfunctions of 
He"*" within a box with free-boundary conditions: eigen¬ 
functions of 

i?+ = -Jv?--+Lb , (3) 

2 ri 

where Lf, is a Bloch operator m- The outer electron is 
described using an FD representation. The presence of a 
Bloch operator for the inner electron leads to boundary 
derivative terms arising at the boundary between the one- 
and two-electron outer regions. As a consequence, the 
TDSE for the one electron outer region becomes 


i^fp{r2,t) = Hfp{r2,t) + 


dt" 


dri 


( 4 ) 


where fp(r 2 ,t) is an FD representation of the wavefunc¬ 
tion at radial distance r 2 and time t in channel p. flpq 
is the surface amplitude that links channel p to the two 
electron finite difference grid r 2 ), and A is the an- 

tisymmetrisation operator. This equation is similar to 
the propagation equation for inner region single ionisa¬ 
tion RMT theory mm- However, the final term on 
the right-hand side connects the wavefunction fp(r 2 ,t) 
in the one-electron outer region with the wavefunction 
^'g(ri,r 2 ) in the two-electron outer region. 

In the inner region, the propagation equations are iden¬ 
tical to those in RMT theory for single ionisation mm- 

= -iy^Hkk'Ck'jt)-^y^ujpk ^ 

k' p r—a 

( 5 ) 

In this equation, the coefficients Ck{t) are coefficients 
of the two-electron R-matrix basis functions in the in¬ 
ner region, and copk are boundary amplitudes of the two- 
electron functions at the boundary r 2 = a. 

RMT theory for single ionisation has previously 
adopted Arnold! propagators for time propagation, re¬ 
quiring separate propagators for the homogeneous TDSE 
and for each boundary term that contributed to the prop¬ 
agation. The complexity arising from multiple propaga¬ 
tors is avoided here through the use of a simpler Taylor 
series computational scheme. 


At each stage of the calculation, we consider the anti¬ 
symmetry of the wavefunction under particle exchange. 
We have implemented this anti-symmetry by considering 
the two-electron outer region wavefunction for all ri > 
a, r 2 > a. The link between the one-electron outer region 
and two-electron outer region is taken along both bound¬ 
aries of the two-electron outer region, ri = a, r 2 > a and 
ri > a, r 2 = a. A phase change can be included to ac¬ 
count for spin symmetry of the electron pair. 

Analysis of the final wavefunction occurs through pro¬ 
jection of the final state onto uncorrelated products 
of eigenfunctions of the one-electron He"*" hamiltonian: 
Fi^{Ei,ri)Fi^{E 2 ,r 2 )- As the two electrons travel far 
from the core and from each other, the energies can be 
interpreted as momenta p squared: = pi • pi/2m = 

kflA/2m, or in a.u.: ki = ^/2Ei. In this communication 
we use k rather than E and generally refer to it as the 
momentum of the electron. 

This approach is correct in the limit in which the 
two ejected electrons are sufficiently weakly interacting. 
In the present calculations we propagate the electrons 
for (typically) 70 field periods after the laser pulse has 
ramped to zero, during which time their Coulomb repul¬ 
sion enhances their spatial separation. We note that k is 
derived from energy and includes contributions from the 
Coulomb potential, but these contributions appear to be 
negligible. If we allow the two-electron wave packet to 
depart the core for durations longer than the usual 70 
field periods, then the calculated energy spectra do not 
change. 

To separate the single-photon yield from the two- 
photon process, we define a cut-off c such that all wave- 
function population of momenta fci, k 2 where + fc| < (? 
is considered to be resulting from the single-photon pro¬ 
cess. Here c is chosen for an individual photon energy E 
as c >> \/2{E — Ip), where Ip is the ionisation poten¬ 
tial for He. From this we find a single-photon double¬ 
ionisation yield a from which we can calculate a corre¬ 
sponding cross section cr^+ as = au! I{t)dt {I{t) 
is the intensity of the pulse at time t and to is the fre¬ 
quency of the pulse). 

The presence of excited bound states (both single¬ 
electron and double-electron) in the final state wavefunc¬ 
tion complicates the analysis of the final-state energy dis¬ 
tributions of the two ejected electrons. As we are cur¬ 
rently only interested in double-ionisation wave packets, 
a Gaussian mask was applied to the wavefunction in the 
regions ri < 35ao, ^2 < 35ao to hide the singly bound 
states and rf-l-r^ < (150ao)^ for the doubly-bound states. 
Care was taken to ensure that by the end of the time- 
propagation, the final yield had stabilised with respect to 
time, indicating that all double-ionisation wave packets 
had propagated into the unmasked area. 

To demonstrate the capabilities of this method, five 
sets of photo-ionisation data were calculated for pulses 
of photon energies 84 eV, 99 eV, 125 eV, 150 eV, and 
180 eV and of peak intensity 4 x lO^^W/cm^. The pulse 
is comprised of 15 cycles (two cycles of sin^ turn-on, 11 
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FIG. 1. (Colour online) Two-electron helium wavefunction 
density (ao~^) 70 cycles after the end of a 15 cycle pulse. 
ri corresponds to the radial distance of electron 1 and r 2 
corresponds to the radial distance of electron 2. The pulse 
has a photon energy of 150 eV and a peak intensity of 10^'* W 
cm“^. All distances are given in bohr radii (ao). Boundaries 
at ri,r 2 = 25ao divide the three RMT regions. 


FIG. 2. (Colour online) Two-electron probability distribution 
in momentum space calculated from the final state wavefunc¬ 
tion (shown in figure[^. fci corresponds to the radial momen¬ 
tum of electron 1 and k 2 corresponds to the radial momentum 
of electron 2. The pulse has a photon energy of 150 eV and 
a peak intensity of 10^“^ W cm~^. All momenta are given in 
a.u.. 


cycles of constant peak intensity, and two cycles of sin^ 
turn-off). A grid spacing of Ax = 0.25ao is used for both 
dimensions of the FD grid, and the boundary between 
the regions is placed at 6 = 25ao. The propagator is a 
6th order Taylor series propagator with a time step of 
0.028as. A basis of 50 B-splines is used to describe the 
single electron functions in the inner region from which 
the two electron states were constructed. Two electron 
states with energies above 1000 a.u. were excluded from 
the calculations. Individual electrons are limited to a 
maximum angular momentum £ = 3 and the atom to 
a maximum angular momentum of L = 2. The two- 
electron outer region is limited to ri + r 2 ^ 900ao. 

Figure [2 shows a two electron wavefunction 70 cycles 
after the end of a 15 cycle laser pulse of photon energy 150 
eV with the momentum transform for the large t limit cal¬ 
culated from this wavefunction given in figurej^ A single¬ 
photon non-sequential process (where both electrons are 
simultaneously ionised by a single photon) is indicated in 
the momentum plot by the arc at {kf + fc^) ~ (2.3a.u.)^, 
as would be predicted given the necessary energy sharing 
for this process. This process is found in the correspond¬ 
ing wavefunction density (figure]^ in the arc in the region 
^ 1+^2 = (220ao)^. In addition, evidence of a two-photon 
sequential process is seen in figure [l] with the outer and 
inner electrons at a distance of « 300ao and « 250ao 
from the nucleus respectively. This process may also be 
seen in the momentum transform for ki and k 2 ~ 3a.u.. 

Processes other than double ionisation are also cal¬ 
culated using this approach, and their effects are visi¬ 
ble in figure Single ionisation is visible close to axes 
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FIG. 3. (Golour online) Single-photon double ionisation cross 
sections calculated using RMT alongside experimental data 
from [22] and theoretical data from |23 | . 


fi,f 2 = 0, and excitation to doubly-excited states is vis¬ 
ible near the nucleus, in addition to the remaining pop¬ 
ulation in the ground state. While further information 
about these processes can in theory be extracted from the 
final wavefunction data, none of these processes appear 
in the momentum transform due to the Gaussian mask 
acting on all non double-ionisation processes. 

Figure shows single photon double-ionisation cross 
sections obtained from a yield calculated using the mo- 

















4 


qQ. 


t A 



FIG. 4. (Colour online) Two photon momentum density 
across an arc corresponding to C^in < fci + fcl < c^ax where 
k\ is the radial momentum of electron 1 and k 2 is the radial 
momentum of electron 2. Arrows indicating the expected an¬ 
gle of peaks caused by a sequential process involving an exci¬ 
tation of the un-ionised electron to 2s or 2p are shown for the 
84 eV and 99 eV data. 


mentum transform of the final wavefunctions. This data 
is given in comparison with experimental data |22| and 
theoretical data [53]. The theoretical RMT data follows 
the overall pattern of the experimental data, with the 
greatest disagreement (« 15%) seen at 84 eV, reducing 
to « 4% at 125 eV, 150 eV and 180 eV. To examine 
the effect of pulse length, cross-sections were calculated 
at photon energies of 84 eV and 99 eV with a 30 cycle 
pulse (4 cycles ramp-on and 4 cycles ramp-off). Simi¬ 
lar cross-sections were calculated for both pulse lengths, 
indicating that the cause of the greater error at these 
photon energies lies elsewhere. 

At low momenta, the double ionised He^"*" wave packet 
is difficult to distinguish from the long tails of higher 
energy He’*' bound states. The calculation of the final 
double ionisation yield for these near-threshold photon 
energies most likely contains a contribution from these 
bound states. Since it is difficult to measure this contri¬ 
bution exactly without propagating over a prohibitively 
large configuration space, we consider this effect to be a 
probable source of error. In addition, the differences be¬ 
tween the current and the benchmark cross-sections may 
be reduced by extending the range of angular momenta 
over which the wave function is represented. 

To demonstrate that the RMT approach has the capac¬ 
ity to extract information about the two-photon sequen¬ 
tial process (as well as the single photon nonsequential 
process), a mask was applied to the momentum trans¬ 
form so that only the momentum density in the region 
Cmm < + k 2 < c^ax the of angular momentum 

couplings L = 0 and L = 2 was retained. Cmin and 
Cmax are given values according to the photon energy 


for the two-photon sequential double ionisation momenta 
(for example for a photon energy of 150 eV in Figure 
Cmin = 3.9a.u. and Cmax = 4.2a.u.). The momentum 
density of the remaining arc is plotted against the an¬ 
gle 9 = tan~^{ki/k 2 ) in figure |4 These spectra are re¬ 
lated to ejected-electron spectra, originally explored for 
sequential double ionisation by [24] . who showed their 
dependence on pulse length. 

The two largest peaks in each photon energy are di¬ 
rectly in the region expected for the sequential process 
where an electron is excited from the ground state into 
the continuum, leaving the bound electron in the He^ Is 
state which is ionised by a later photon. In the data for 
84 eV and 99 eV features an order of magnitude smaller 
are visible, which occur at momenta corresponding to a 
two photon process (shown by the arrows in figure 
where the bound He'*' electron is excited to either the 
2s or 2p state before being ionised by the second photon 
(as discussed in detail in |25l|26]). While these processes 
should also be present in the 125 eV, 150 eV and 180 
eV spectra, they occur at angles where the Is process 
dominates, making them difficult to observe. The 30 cy¬ 
cle pulse data shows these features more distinctly than 
the 15 cycle pulse. For these longer pulses, an additional 
minimum is seen in the 84 eV spectra, corresponding to 
a side-band caused by the pulse length. 

In conclusion, we have combined a two-electron FD 
outer-region approach with the RMT method. This al¬ 
lows the double-ionised wave-packet to be propagated 
over a larger configuration space. The accuracy of the 
approach is demonstrated by the determination of Ffe 
single-photon double ionisation cross sections for photon 
energies in the region 80 eV to 180 eV. We obtain agree¬ 
ment with experiment and existing theory to within 15% 
near the single photon double ionisation threshold, and 
to within 4% for higher photon energies. The capabil¬ 
ity to investigate higher-order processes is demonstrated 
through the observation of signatures associated with a 
sequential ionisation process involving excited states of 
the intermediate He’*' ion. 

This agreement demonstrates the feasibility of attach¬ 
ing a two-electron FD region to an R-Matrix outer region 
using the RMT methods, and that this approach can be 
applied to predict a wide range of experimental observ¬ 
ables. The FD method is highly parallelisable, and the 
current program scales linearly with the area covered by 
the two-electron double outer region. 

The RMT method has the potential to study double¬ 
ionisation with full-correlation in general atoms. In order 
for this potential to be realised, it will be necessary to 
develop a multi-electron inner region basis set with full 
correlation to model double ionisation in the inner region, 
in addition to a corresponding outer region basis set. 
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